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基于 FLUS 模型 的 天 山 山 区 未 来 土地 利用 变化 预 估 


Ram”, 


摘 要 : 基于 天 山 山 区 53 个 气象 站 点 ,31 个 CMIP5 全 球 气 候 模式 ,集合 逐 入 
等 资料 ,利用 FLUS 模型 , 预 估 了 2050 年 RCP2.6、RCP4.5 和 RCP8.5 排放 情景 下 天 山 山区 土地 利用 变化 。 
明 :Q@ FLUS 模型 能 够 较 好 地 模拟 研究 区 历史 土地 利用 变化 ,具有 较 高 的 模拟 精度 。@ 与 基准 
相 比 ,3 种 情景 下 多 年 平均 气温 升 高 了 1.67 ~2.16 Y ,其 


多 年 平均 降水 量 增加 22 ~25 mm ,3 种 情景 下 个 别 


ae 441,2 
李雪梅 “， 
(1. 兰州 交通 大 学 测绘 与 地 理 信息 学 院 ,甘肃 兰州 730070; 2. 甘肃 省 地 理 国 情 监测 工程 实验 室 ,甘肃 兰州 

3. 济南 市 交通 局 运输 管理 办 公 室 ,山东 济南 


陈 庆 伟 ， 田 亚 林 “ 


730070 ; 
250014 ) 


TF 


平均 气候 资料 ,社会 经 济 和 自然 地 形 
结果 表 
期 (1970 一 1999 年 ) 


P 春 季 和 夏季 增 温 幅 度 较 小 ,秋季 和 冬季 增 温 幅 度 较 大 。 
年 份 较 基准 期 出 现 降水 减少 现象 ,但 总 体 呈 增加 趋势 。@) 2050 
年 RCP2.6、RCP4.5 和 RCP8.5 的 预 佑 结果 与 2017 年 实测 土地 利 
34.6% 38.4% 。 前 两 种 情景 城市 用 地 面积 增幅 不 大 ,部 分 耕地 转换 为 城市 


性 结 果 相 比 ,天 山 山区 冰雪 面积 分 别 减 少 29.2% 、 
用 地 和 草地 ,林地 面积 少量 增加 。 而 


RCP8.5 情景 下 城市 用 地 面积 为 2017 年 的 近 3 倍 , 部 分 林地 转换 为 草地 。3 种 情景 下 水 体 总 面积 略 有 增加 ,未 利用 


地 时 减少 趋势 。 
关键 词 : 气候 变化 ; 情景 模拟 ; 土地 利 


在 全 球 变 暖 的 背景 下 ,气候 变化 对 全 球 的 土地 
利用 类 型 .生态 环境 和 社会 经 济 产 生 了 一 系列 影 
me? 。 气 候 变 化 (全 球 变 暖 ,极端 天 气 等 ) 和 生态 退 
化 (水 文 变化 .土壤 侵蚀 等 ) 对 自然 景观 的 动态 变化 
具有 长 期 的 影响 ” 。 而 土地 是 人 类 赖 以 生存 的 基 
础 ,土地 利用 /覆盖 变化 (land-use and land-cover 
change, LUCC) 对 自然 环境 和 人 类 生存 发 展 有 着 深 
远 的 影响 。 研 究 LUCC 机 理 、 预 测 土地 利用 演变 趋 
势 , 有 利于 合理 开发 利用 土地 资源 和 实现 土地 利用 
的 可 持续 发 展 ” 。 

利用 全 球 气 候 模 式 (global climate model， 
GCM ) 对 未 来 气候 的 情景 预 估 成 为 世界 研究 的 热 
点 。2013 年 ,联合 国政 府 间 气 候 变化 专门 委员 会 
(Intergovernmental Panel on Climate Change, IPCC) 
发 布 了 第 5 次 评估 报告 ,该 报告 中 未 来 气候 变化 
的 预 估 主 要 基于 全 球 耦 合 模型 比较 计划 第 五 阶段 
(Coupled Model Intercomparison Project Phase 5, 
CMIPS ) 的 46 个 地 球 系统 模型 结果 ,介绍 了 情景 和 
模型 的 不 确定 性 并 模拟 了 21 世纪 及 以 后 气候 变化 
的 预测 结果 中 。CMIP5 模式 数据 采用 了 新 的 温室 
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变化 ; FLUS 模型 ; 天 山 山区 


气体 排放 情景 , 称 为 "典型 浓度 目标 ”(representative 
concentration pathways ，RCPs ) EO 。 为 了 消除 单 
一 模式 的 偏差 ,采用 多 模式 集合 方法 得 到 的 结果 和 
观测 结果 最 相似 ,更 接近 于 真实 的 自然 状态 。 
影响 LUCC 的 因素 众多 ,致使 土地 利用 类 型 转 
换 过 程 极其 复杂 ,构建 土地 利用 模型 有 助 于 对 复杂 
的 系统 进行 简化 ,使 人 们 能 够 更 加 清晰 地 理解 土地 
利用 变化 的 过 程 , 并 对 未 来 的 变化 做 出 预 估 ”" o 
由 Liu 等 中 提出 的 FLUS(future land use simulation ) 
模型 ,能 够 较 好 地 处 理 各 驱动 力 因子 之 间 的 关系 , 模 
拟 各 土地 利用 类 型 之 间 的 相互 作用 ,具有 和 较 高 的 模 
拟 精度 ,对 比 了 FLUS 模型 与 其 他 土地 利用 模拟 模 
型 的 优 劣 ,并 模拟 了 IPCC AR4 四 种 气候 变化 情景 
下 中 国土 地 利用 动态 变化 ,分 析 了 影响 土地 利用 类 
型 转变 的 因子 ; 周 浩 等 "通过 遥感 影像 分 析 得 出 ， 
1990—2014 年 三 江平 原 典 型 流域 旱地 面积 变化 幅 
度 较 小 ,水 田 面积 持续 增加 ,并 运用 FLUS 模型 预 估 
了 RCPs 情景 下 耕地 动态 变化 特征 ; 张 子 明 等 "” 利 
用 FLUS 模型 对 抚州 市 临 川 区 的 城镇 用 地 增长 边界 
划 定 进行 了 研究 ,并 对 2025 年 研究 区 城镇 用 地 变化 
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进行 了 情景 模拟 。 

近 50 a 来 (1961 一 2010 年 ) 西 北 干旱 半 干 旱地 
区 平均 增 温 幅 度 达 1. 65 Cm ,新疆 天 山 山区 被 誉 
为 “中 亚 水 塔 " 和 “固体 水 库 ”, 由 于 其 独特 的 山 ` 盆 、 
谷地 相间 的 地 形 地 貌 特 征 , 该 区 域 对 全 球 气 候 变 化 
响应 尤为 敏感 ""。 鉴 于 此 ,本 文选 取 天 山 山区 作为 
人 研究 区 ,其 于 FLUS 模型 ,结合 CMIPS 不 同 RCP HE 
放 情 景 下 的 31 个 全 球 气候 模式 ,集合 年 平均 气候 资 
BL .社会 经 济 和 自然 地 形 等 资料 , 预 估 了 天 山 山 区 
2050 年 的 土地 利用 变化 。 本 研究 可 为 天 山 山 区 未 
来 生态 环境 管理 .土地 资源 合理 配置 以 及 洪涝 灾害 
防治 提供 科学 的 决策 依据 。 


1 资料 和 方法 


1.1 研究 区 概况 

我 国境 内 的 天 山 山区 横贯 新 疆 中 部 , 占 新 疆 总 
面积 的 1/3 ,其 北面 为 准噶尔 盆地 ,南面 为 塔里木 盆 
地 。 天 山 将 新 疆 分 为 北 疆 和 南 疆 ,是 南北 疆 气 候 的 
分 水 岭 !。 东 西 绵延 1 700 km ,南北 宽 250 ~ 350 
km, 总 面积 约 为 5.7 x10” km ,山地 的 平均 海拔 高 
度 为 4000 m, 

天 山 山区 深 居 我 国 西北 内 陆 ,属于 干旱 区 ,具有 
典型 的 大 陆 性 气候 特征 。 由 于 降水 受 天 山 独 特 的 地 
形 地 貌 和 西风 环流 的 影响 ,东部 地 区 少 于 西部 , 南 坡 
少 于 北 坡 ,平原 少 于 山区 ,背风 坡 少 于 迎风 坡 ; 降 水 
随 着 海拔 的 升 高 先 增多 后 减少 ,最 大 降水 带 一 般 在 
中 高 山地 区 ,盆地 为 少 十 区。 而 天 山 山 区 的 年 均 气 
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温 西部 地 区 低 于 东部 , 同 海拔 北 坡 低 于 南 坡 。 受 海 
影响 天 山 具 有 明显 的 气候 垂直 地 带 性 ,形成 了 不 

同 植被 和 土壤 的 自然 景观 " 。 

虽然 天 山 位 于 炎热 干燥 的 干旱 区 ,高 大 的 山地 
却 终 年 寒 凉 湿润 , 积 雪 、 冰 川 与 冻 土 分 布 广 泛 , 冰 和 雪 
的 消融 在 一 定 程度 上 改变 了 天 山南 北 某 些 地 理 要 素 
的 状况 , 极 大 地 影响 了 该 地 区 人 民 的 生产 和 生 
wg”, 
1.2 数据 资料 

本 文 所 用 数据 包括 天 山 山 区 53 个 气象 站 点 
(图 1)1970 一 1999 年 的 年 平均 气温 和 年 累积 降水 
资料 ;31 个 CMIPS 全 球 气 候 模式 集合 的 RCP2. 6, 
RCP4.5 和 RCP8.5 排放 情景 下 2021—2050 年 的 年 
平均 气温 和 年 累积 降水 资料 (http :// www. iiasa. ac. 
at/web-apps/tnt/RepDb/dsd? Action = htmlpage& 
page = welcome ) , CMIPS 全 球 气 候 模式 的 详细 介绍 
参见 https://pemdi. nl. gov/。 男 外 ,所 用 数据 还 包 
括 2001—2017 年 空间 分 辨 率 为 500 m FY) EHA ate 
类 型 MCD12Q1 产品 数据 ,该 数据 来 源 于 美国 航空 
航天 局 (NASA https://urs. earthdata. nasa. gov/) ; 
全 国 基础 地 理 信息 数据 来 源 于 全 国 地 理 信息 资源 目 
SEARS KB (http://www. webmap. cn/ ) ;1 km 分 辨 
率 人 口 和 GDP 数据 ,社会 经 济 数据 来 源 于 中 国 科学 
院 资源 环境 科学 数据 中 心 (http://www. resde. 
cn/) 。 
1.3 温室 气体 排放 情景 

RCP2. 6 是 温室 气体 排放 和 辐射 强迫 最 低 的 情 


300 km 


1 研究 区 气象 站 点 空间 分 布 


Fig.1 Distribution of the meteorological stations in the study area 


于 g 


景 ,辐射 强迫 先 增加 后 减少 ,至 2100 年 减少 为 2.6 
W- mm ,全 球 平均 温度 上 升 限 制 在 2 % 以 内 。 
2010 一 2100 年 累计 温室 气体 排放 比 基 准 期 减少 
70% ,GDP、 人 口 增长 率 和 技术 创新 为 最 低 水 平 , 促 
进 生物 质 能 利用 ,提倡 恢复 森林 。 

RCP4.5 为 中 等 排放 情景 ,至 2100 年 辐射 强迫 
为 4.5 W. m“。 该 情景 考虑 了 与 全 球 经 济 发 展 相 
适应 的 温室 气体 等 其 他 物质 的 排放 ,用 最 低 代 价 达 
到 辐射 强迫 目标 。 技 术 发 展 较 快 ,通过 多 用 电能 \ 低 
排放 能 源 技术 限制 温室 气体 的 排放 ,经 济 .社会 和 环 
境 可 持续 发 展 ' 。 

RCP8.5 是 最 高 的 温室 气体 排放 情景 ,至 2100 
年 辐射 强迫 升 至 8.5 Wem? KARAOKE, 
技术 创新 不 高 能源 改善 缓慢 ,缺少 应 对 气候 变化 的 
政策 ,大 规模 的 人 类 活动 加 速 温室 气体 排放 ,增加 了 
气温 和 降水 的 异 质 性 。 

由 于 大 多 数 模 型 未 能 在 RCP6. 0 情景 进行 预 估 
试验 ,本 文 没 有 对 该 情景 中 的 仿真 结果 进行 预测 ,此 
SE TTA NESE AS BCBS 
1.4 RCPs 情景 GDP 和 人 口 增长 率 

根据 新 疆 统计 年 鉴 . 国 家 统计 局 的 有 关 课 题 组 
对 “十 三 五 ”时 期 的 经 济 增长 率 测算 和 模拟 ,以 及 国 
务 院 新 闻 办 公 室 2000 年 发 表 的 《中 国 21 世纪 人 口 
与 发 展 》 白 皮 书 和 中 国 第 六 次 人 口 普查 数据 ,在 新 
ii 2007—2017 年 的 平均 GDP 和 人 口 增长 率 的 基础 
上 进行 微调 整 , 设 定 RCPs 情景 下 2017 一 2050 年 的 
GDP 和 人 口 增长 率 ( 表 1)。 


表 1 RCPs 情景 GDP 和 人 口 增长 率 
Tab.1 Growth rates of GDP and population under 


RCPs scenario 


驱动 因子 RCPs 情景 增长 率 

GDP RCP2. 6 从 14% 线 性 下 降 到 4% 
RCP4.5 从 14% 线 性 下 降 到 7% 
RCP8.5 保持 14% 增长 

AW RCP2.6 从 6%o 线 性 回落 到 4%o 
RCP4.5 保持 5%o 增 长 
RCP8.5 保持 6%o 增 长 

1.5 全球 气 候 模式 降 尺 度 


CMIP5 全 球 气候 模式 模拟 精度 依旧 较 低 ,忽略 
了 小 太 度 上 由 于 地 表 形 态 .生物 变化 和 地 表 物 理 过 
程 等 因素 产生 的 差异 。 一 般 通 过 降 尺 度 的 方法 
来 降低 误差 ,获取 更 高 精度 的 气候 资料 。 本 文通 过 


双 线 性 插值 法 将 31 个 气候 模式 1° x 1? 的 年 平均 数 
据 插 值 到 站 点 ,以 1970 一 1999 年 为 基准 期 ,通过 偏 
差 修正 中 的 变化 因子 法 (Delta ) 对 未 来 2021—2050 
年 各 情景 逐年 数据 进行 了 降 尺 度 。 

Delta 方法 是 一 种 比较 简单 但 在 GCMs 数据 降 
尺度 时 常用 的 方法 。 降 水 选择 的 是 相对 变化 量 , 即 
比较 GCMs 输出 格 网 未 来 各 情景 模拟 某 年 累积 降水 
量 与 未 来 各 情景 模拟 多 年 平均 累积 降水 量 ,计算 各 
网 格 降水 量 的 变化 率 ,各 个 气象 站 基准 期 实测 多 年 
平均 累积 降水 量 乘 以 降水 量 的 变化 率 ,得 到 未 来 各 
情景 年 累积 降水 量 数据 。 气 温 选 择 的 是 绝对 变化 
量 ,计算 GCMs 输出 格 网 未 来 不 同情 景 某 年 平均 气 
温 与 未 来 多 年 平均 气温 的 变化 量 ,将 该 变化 量 加 到 
网 格 内 各 个 气象 站 基准 期 实测 多 年 平均 气温 上 ,得 


到 未 来 各 情景 年 平均 气温 数据 。 
Delta 方法 处 理 降 水 数据 计算 公式 为 : 
Por 
neh (1) 
Delta 方法 处 理气 温 数据 计算 公式 为 : 
T,=T,+(Ta- Te) (2) 


式 中 :P, 为 实测 多 年 平均 累积 降水 量 ; Pe 为 模拟 未 
来 多 年 平均 累积 降水 量 ;Pe 为 模拟 未 来 某 年 累积 降 
水 量 ;P 为 Delta 方法 得 到 的 未 来 年 累积 降水 量 ; 7， 
为 实测 多 年 平均 气温 ; Te 为 模拟 未 来 多 年 平均 气 
温 ; Tc 为 模拟 未 来 某 年 平均 气温 ;Tr 为 Delta 方法 得 
到 的 未 来 年 平均 气温 。 

1.6 数据 处 理 

FLUS 模型 中 的 数据 资料 必须 转换 成 相同 的 分 
辨 率 和 投影 。 本 文 将 所 有 数据 分 辨 率 重 采样 为 
1 km ,投影 坐标 系 为 WCGCS_1984_UTM_zone_ 43N。 

在 ArcGIS 中 ,将 31 个 气候 模式 年 均 气 温和 年 
累积 降水 的 降 尺 度数 据 , 采 用 反 距 离 权 重 法 插值 为 
1 km 机 格 图 像 ;将 30 m 的 DEM 数据 , 重 采样 为 
1 km 空间 分 辨 率 , 再 提取 坡度 和 坡 向 ;利用 欧 氏 距 
离 分 析 全 国 基础 地 理 信息 数据 中 的 镇 点 .道路 、 铁 
路 河流 数据 ,得 到 km 距离 图 (图 2) ;采用 栅 格 计 
算得 到 RCPs 情景 下 GDP 和 人 口 图 。 

MODIS 土地 利用 数据 MCD12Q1 采用 IGBP( Tn- 
ternational Geosphere-Biosphere Programme ) 的 全 球 植 
被 分 类 方案 ,本 文 在 ArcGIS 中 通过 重 分 类 分 为 水 、 
林地 草地 耕地 ,城市 ` 冰 雪 和 未 利用 土地 7 类 ( 表 
2) , 重 采样 为 1 km 空间 分 辩 率 。 
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O) 到 铁路 距离 
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0.14 元 .lm- 
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0m 
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Om — | 


图 2 土地 利用 变化 驱动 力 因 子 


Fig.2 Driving factors of land use change 


表 2 土地 利用 类 型 分 类 
Tab.2 Classification of land use types 
本 文 分 类 ICBP 全 球 植被 分 类 本文 分 类 1GBP 全 球 植被 分 类 


水 水 草地 木 本 热带 稀 树 草原 
永久 湿地 热带 稀 树 草原 
林地 常 绿 针 叶 林 草地 
常 绿 阔 叶 林 耕地 农用 地 
落叶 针 叶 林 农用 地 /自然 植被 
落叶 阔 叶 林 城市 城市 和 建筑 区 
混交 林 冰雪 冰雪 
稠密 灌 从 未 利用 土地 未 利用 土地 
Pi HE DA 


1.7 FLUS 模型 

FLUS 模型 由 2 个 主要 模块 构成 ,分 别 为 基于 神 
经 网 络 (ANN ) 的 适宜 性 概率 计算 模块 和 基于 自 适 
应 惯性 机 制 的 元 胞 自动 机 模块 ,ANN 分 为 3 层 , 第 1 
层 是 输入 层 ,第 2 层 是 隐藏 层 ,第 3 层 是 输出 层 。 首 
先 ,将 一 期 土地 利用 数据 结合 人 类 社会 经 济 与 自然 
效应 的 多 种 驱动 力 因 子 (气温 、 降 水 、 地 形 、 人 口 分 
布 等 ) 中 采用 ANN 整合 并 计算 研究 区 域内 每 种 土 
地 利用 类 型 在 每 个 像 元 上 的 出 现 概 率 , 获 取 研 究 区 
各 类 土地 的 适宜 性 概率 。 其 次 ,基于 轮 盘 赌 选择 的 
自 适应 惯性 竞争 机 制 ,在 模拟 过 程 中 可 以 有 效 地 处 


理 人 类 活动 和 自然 因素 影响 下 各 种 土地 利用 类 型 转 
变 的 不 确定 性 和 复杂 性 ,获取 研究 区 土地 利用 模拟 
图 "2 。 其 表达 公式 为 : 

p(k,t,l) = 之 iv x sigmoid[ net; (k,l) | 


1 
= wi, “了 $e HED (3) 


AP p(k, t, DRRR t KEME kT 
EAR w, sigmoid 分 别 为 隐藏 层 与 输出 层 间 
的 权 值 和 激励 函数 ;net; (上 ,7) 表示 在 第 i 个 隐藏 层 


栅 格 在 时 间 1 的 响应 值 。 
Lethe) =1 (4) 
式 中 :适宜 性 概率 p(k,t,7) RIRE LAL EY A k 
上 ,各 类 的 适宜 性 概率 之 和 为 1。 
Inertia! i |M |< |M | 
I M,” 0> M7? sm"! 
ieii = nertiay git 7 of >M, >M, (5) 


1-1 


are 1-2 
Inertia, aTe ka >M,” >0 
t 


式 中 :Inertia 为 土地 需求 数量 与 土地 数量 的 差异 决 


干旱 区 


定 自 适应 惯性 系数 ;MM 分 别 表示 i 种 用 地 类 
型 在 1- 1 7-2 两 个 不 同时 刻 的 机 格 数 与 所 需求 数 
Ze st 种 地 类 在 1 时 刻 的 自 适应 惯性 系数 表示 为 
Inertia’ 。 


TPS, , =p(k,t,l) x Q),, x Inertia, x (1 -sc,,,) (6) 


IUP TPS, dean LAA HA k EEJ t 种 用 地 类 
7 EN) EER ; scn N EMA FRY c 转化 为 类 型 ;的 
成 本 ;1 -sc, ,表示 发 生 转化 的 难 易 程度 ;Q; ,为 领域 
作用 。 


onl cy ! =t) 


L 一 一 一 
Q, = NxN-1 w, (7) 


式 中 : X con(c | =t) RIRE N xN 的 摩尔 邻 域 窗 


口 , 上 一 次 迭代 结束 后 第 i 种 地 类 的 栅 格 总 数 ,本 文 
N =3 3, 为 邻 域 影响 程度 的 权重 。 


2 结果 与 分 析 


2.1 数据 检验 和 模型 验证 

2.1.1 降 尺 度数 据 检验 ”检验 CMIP5 多 模式 集合 
基准 期 的 气温 和 降水 数据 在 天 山 山区 降 尺 度 的 结 
见 表 3。CMIP5 多 模式 集合 对 于 天 山 山区 基准 期 的 
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气温 模拟 偏 低 ,降水 量 模拟 偏 高 ,但 是 冬季 气温 和 夏 
季 降 水 模拟 值 与 实测 值 误差 较 小 。 经 过 Delta BER 
度 处 理 后 基准 期 的 多 年 平均 值 与 实测 值 相近 ,其 中 
气温 ( 除 冬 季 ) 和 降水 较 实测 值 略微 偏 低 ,气温 误差 
为 0.15 左右 ,降水 相对 误差 都 在 1% ~5% 。 

2.1.2 FLUS 模型 精度 检验 ”采用 2001 年 土地 利 
用 现状 图 模拟 2010 年 土地 利用 对 FLUS 模型 进行 
率 定 , 并 用 2010 年 土地 利用 现状 图 模拟 2017 年 土 
地 利用 进行 验证 。 首 先 ,将 2001 年 土地 利用 现状 图 
结合 2001—2010 年 驱动 力 因 子 作 为 输入 层 , 输 出 层 
有 7 个 神经 元 (对 应 7 种 土地 利用 类 型 ) ,每 种 土地 
利用 类 型 均匀 取样 2% 进行 神经 网 络 训练 ,获取 各 
类 土地 的 适宜 性 概率 图 。 然 后 ,将 2001 年 土地 利用 
图 和 土地 利用 类 型 适宜 性 概率 输入 自 适应 惯性 机 制 
的 元 胞 自动 机 模块 ,运行 得 到 2010 年 土地 利用 模拟 
图 。 最 后 ,将 2010 年 土地 利用 现状 图 与 模拟 图 进行 
对 比分 析 ( 图 3) ,计算 结果 Kappa 系数 为 0.71, 整 
体 精度 为 0.83。 同 上 操作 ,采用 2010 年 土地 利用 
现状 图 模拟 2017 年 土地 利用 ,并 与 2017 年 土地 利 
用 现状 图 对 比 , 计 算 结果 Kappa 系数 为 0.74 ,整体 
精度 为 0.84。 由 此 表明 ,FLUS 模型 能 较 好 地 模拟 
出 土地 利用 类 型 分 布 ,该 模型 可 用 于 土地 利用 变化 


表 3 气温 和 降水 降 尺度 前 后 对 比 结果 


Tab.3 Compared results of temperature and precipitation before and after downscaling 


年 份 ”观测 /CMIP5 模式 预 估 ates a eerie Wap 
年 £ E K 3% 年 春 夏 秋 冬 
1970 一 1999 观测 7.22 9.41 21.50 7.49 -9.52 184. 56 198. 48 326. 52 147.48 65.76 
降 尺 度 前 4.63 5.46 19. 36 4.24 一 10.S2 261.12 347.88 344. 28 193.56 158.76 
降 尺 度 后 7.09 9.23 21.11 7.36 一 9.34 181. 20 194. 88 320.76 144. 84 68.76 
N (a) 200148 (b) 2010 年 (c) 2017 年 


(d) 2010 年 


(e) 2017 年 


E 
[未 利用 土地 


图 3 土地 利用 现状 与 FLUS 模拟 对 比 


Fig.3 Comparison between present situation of land use and FLUS simulation 
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的 预 估 。 
2.2 气温 和 降水 的 变化 预 估 

CMIPS 多 模式 集合 在 RCP2. 6、RCP4. 5 和 

RCP8.5 情景 下 ,天 山 山 区 的 年 平均 气温 和 年 累积 
降水 量 经 过 Delta 方法 处 理 后 ,以 1970 一 1999 年 作 
为 基准 期 , 预 估 2021—2050 年 相对 于 基准 期 的 变 
化 。 
2.2.1 年 平均 气温 变化 分 析 3 种 情景 下 ,天 山 山 
区 2021—2050 年 CMIP5 多 模式 集合 多 年 平均 气温 
较 基 准 期 均 有 升 高 ,RCP2. 6, RCP4. 5 和 RCP8.5 情 
景 分 别 增加 了 1.67 °C .1.77% 和 2.16%( 图 4)。 
从 季节 来 看 ,春季 和 夏季 升温 幅度 较 小 ,秋季 和 冬季 
升温 幅度 较 大 。 

RCP2. 6 情景 下 2030 年 .2040 年 .2050 年 平均 
气温 较 基准 期 分 别 增加 1.28 °C .1.26% .0.81 C; 
RCP4.5 情景 下 分 别 增加 1.11 °C 1.29 %C 1.39 C; 
RCP8.5 情景 下 分 别 增加 1.39 % 1.52 € 1.97 C, 
2.2.2 年 累积 降水 变化 分 析 3 种 情景 下 ,天 山 山 
区 2021 一 2050 年 CMIP5 多 模式 多 年 平均 降水 量 较 
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图 4 RCPS 较 基准 期 年 平均 气温 变化 
Fig.4 Amplitude of mean annual temperature during 
2021 — 2050 under RCPs scenario compared with 
that during base period (1971 - 1999 ) 


N (a) RCP2.6 情 景 


(b) RCP4.5 情 景 


基准 期 均 有 增加 。 由 于 天 山 山 区 地 处 干旱 区 ,整体 
增幅 不 大 ,为 22 ~25 mm。 其 中 RCP2.6 情景 增幅 
最 大 ,RCP4.5 情景 增幅 最 小 。 从 季节 上 来 看 ,夏季 
降水 增幅 较 大 ,冬季 降水 量 增 幅 较 少 。3 种 情景 下 
个 别 年 份 较 基准 期 出 现 降水 减少 的 现象 ,但 总 体 呈 
增加 趋势 (图 5)。 

RCP2. 6 情景 下 2030 年 .2040 年 .2050 年 累积 
降水 量 较 基准 期 分 别 增 加 40.9 mm .26.7 mm .90.7 
mm;RCP4.5 情景 下 分 别 增加 37.2 mm.、23.8 mm, 
74.6 mm; RCP8. 5 情景 下 分 别 增加 30. 1 mm 、26.6 
mm 88 mm, 

2.3 ”未 来 土地 利用 变化 分 析 

将 2017 年 土地 利用 现状 图 和 3 种 情景 下 的 驱 
动力 因子 作为 FLUS 模型 的 输入 ,运行 模型 后 输出 1 
km 空间 分 辩 率 的 土地 利用 栅 格 图 ,分析 2050 年 
RCPs 情景 下 天 山 山区 土地 利用 变化 (图 6)。 

2.3.1 RCP2.6 情 景 下 未 来 土地 利用 变化 RCP 
2.6 情 景 是 相对 保守 的 发 展 策略 ,与 2017 年 相 比 ， 
2050 年 水 域 .林地 、 草 地 和 城市 面积 分 别 增加 了 
8.6% 、3% 、10.1% 和 24.7% , 而 耕地 、 冰 雪 和 未 利 


100 p RCP2.6 =m CP4.5 mmm ROP8.5 
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图 5 ROPs 较 基准 期 年 累积 降水 量变 化 
Fig.5 Amplitude of cumulated annual precipitation during 
2021 -2050 under RCPs scenario compared with that 
during base period (1971 - 1999) 


(c) RCP8.5 情 景 


Bit Es Et Exes 草地 [| 未 利用 土地 Eee 
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6 RCPs 情景 下 2050 年 土地 利用 变化 
Fig.6 Land use change in 2050 under RCPs scenario 


用 土地 面积 分 别 减 少 了 24.3% 29.2% 和 3.1%( 表 
4)。 在 不 同 土地 利用 类 型 向 其 他 地 类 转 入 的 过 程 
中 ,耕地 和 未 利用 地 向 草地 转换 的 面积 较 大 ,分 别 为 
13 339 km 和 8 501 km , 占 草 地 转 入 总 面积 的 
57.8% 和 36.8% ;冰雪 转 入 草地 面积 为 1 034 km’, 
占 冰 雪 总 面积 的 24. 6% 。 草 地 和 未 利用 地 土地 利 
用 变化 活跃 ,向 各 类 土地 利用 类 型 均 有 转换 。 城 镇 
周围 耕地 大 面积 地 向 草地 转换 , 奇 台 一 木 件 一 吉木 
BEAR .库尔勒 一 在 普 和 哈密 等 地 区 耕地 明显 减少 ; 城 
市 的 扩张 也 侵占 了 部 分 耕地 ,如 :乌鲁木齐 .吐鲁番 
和 托 克 逊 等 ;其 他 类 型 的 土地 转换 为 草地 的 面积 小 
且 较 为 分 散 。 由 此 可 见 ,RCP2.6 情景 低 人 口 .GDP 
增长 率 和 提倡 恢复 森林 的 政策 ,使 得 天 山 山区 粮食 
压力 减少 ,导致 耕地 面积 大 幅 减少 , 而 林地 略 有 增 
加 。 气 候 变 暖 导致 雪线 上 升 ,2050 年 冰雪 面积 将 减 
少 29.2% ,消融 的 冰雪 转换 为 草地 和 未 利用 地 。 低 
海拔 地 区 范围 较 小 的 冰雪 将 会 消失 ,高 海拔 地 区 冰 
雪 面 积 消融 较 少 。 降 水 的 增多 和 冰雪 消融 的 补给 ， 
导致 艾 比 湖 和 博 斯 腾 湖 面积 有 所 增加 ,水 体 总 面积 


研 究 
增加 。 
2.3.2 RCP4.5 情景 下 未 来 土地 利用 变化 在 


RCP4.5 情景 下 ,2050 年 较 2017 年 水 域 . 林 地 、 草 地 
和 城市 面积 分 别 增加 了 2. 1% 17% . 16. 3% 和 
36.8% ,而 耕地 冰雪 和 未 利用 土地 面积 分 别 减少 了 
29.8% .34.6% 和 7.1% ( 表 5) 。 在 不 同 土地 利用 类 
型 向 其 他 地 类 转 入 的 过 程 中 ,耕地 和 未 利用 地 向 草 
地 转换 的 面积 较 大 ,分 别 为 12 942 km? 和 24 821 
km? , 占 草地 转 和 人 总 面积 的 33.7% 和 64.6% ;草地 转 
换 为 林地 的 面积 为 339 km’, 占 林 地 转 和 人 总 面积 的 
86.3% ;草地 转换 为 耕地 的 面积 为 2 981 km , 占 耕 
地 转 入 总 面积 的 94.7% 。 气 温 的 上 升 和 降水 的 增 
加 ,导致 草地 和 未 利用 地 变化 活跃 ,向 各 类 土地 利用 
类 型 均 有 转换 。 草 地 的 分 布 区 域 集聚 化 ,草地 的 景 
观 破 碎 度 有 所 减少 , 阿 合 奇 西 部 和 巴里 坤 地 区 的 未 
利用 地 转换 为 草地 尤为 明显 ;阿克苏 一 乌 什 、 库 车 和 
库尔勒 等 区 域 的 耕地 明显 减少 并 转换 为 草地 。 城 市 
的 发 展 导 致 部 分 耕地 转换 为 建设 用 地 , 主要 在 乌 鲁 
木 齐 石河子 和 库尔勒 等 地 区 。 可 持续 发 展 的 政策 


表 4 2017 年 与 RCP2.6 情景 2050 年 土地 利用 变化 面积 转移 矩阵 


Tab.4 Transition matrix of land use change in 2017 and 2050 under RCP2.6 scenario /km? 
上 地 利用 类 型 
水 域 林地 草地 耕地 城市 冰雪 未 利用 地 2017 年 合计 
2017 年 水 域 2 815 0 57 1 0 0 23 2 896 
林地 0 1 102 97 2 0 0 2 1 203 
草地 287 93 190 501 1 463 329 2 1 260 193 935 
耕地 16 14 13 339 36 116 318 0 25 49 828 
城市 0 0 49 103 1 841 0 14 2 007 
冰雪 0 24 1 034 0 0 2 846 298 4 202 
未 利用 地 27 6 8 501 41 14 129 220 497 229 215 
2050 年 合计 3 145 1 239 213 578 37 726 2 502 2 977 222 119 483 286 
#5 2017 年 与 RCP4.5 情景 2050 年 土地 利用 变化 面积 转移 矩阵 
Tab.5 Transition matrix of land use change in 2017 and in 2050 under RCP4.5 scenario /km? 
naan RCP4.5 情景 2050 年 
水 域 林地 草地 耕地 城市 冰雪 未 利用 地 2017 年 合计 
2017 年 水 域 2 312 13 408 11 0 6 146 2 896 
林地 3 1 015 164 0 0 9 12 1 203 
草地 199 339 187 097 2 981 105 22 3 192 193 935 
耕地 4 0 12 942 31 813 603 0 4 466 49 828 
城市 0 0 24 27 1 940 0 16 2 007 
冰雪 7 2 67 0 0 2 478 1 648 4 202 
未 利用 地 431 39 24 821 129 98 233 203 464 229 215 
2050 年 合计 2 956 1 408 225 523 34 961 2 746 2 748 212 944 483 286 


秦 圭 瑞 等 :基于 FLUS 模型 的 天 山 山区 未 来 土地 利用 变化 预 估 


使 林地 面积 增加 ,增加 的 区 域 比 较 分 散 。 雪 线 逐 渐 
上 升 ,冰雪 面积 逐渐 减 小 ,更 多 小 范围 低 海拔 的 冰雪 
消失 ,减少 的 冰雪 面积 转换 为 未 利用 土地 和 草地 。 
水 体面 积 变化 同 RCP2. 6 情景 相似 。 

2.3.3 RCP8.5 情景 下 未 来 土地 利用 变化 RCP 
8.5 是 温室 气体 排放 最 高 的 情景 ,高 排放 的 发 展 策略 
导致 很 多 土地 利用 类 型 发 生 了 重大 改变 。2050 年 
较 2017 年 水 域 . 草地、 耕地 和 城市 面积 分 别 增加 了 
5.5% 2% .7.9% 和 199. 6% ,而 林地 冰雪 和 未 利 
用 土地 面积 分 别 减少 了 18.5% .38.4% 和 4.5% (R 
6) 。 在 不 同 土地 利用 类 型 向 其 他 地 类 转 和 人 的 过 程 
中 ,耕地 和 未 利用 地 向 草地 转换 的 面积 较 大 ,分别 为 
2 072 km 和 7 125 km , 占 草地 转 入 总 面积 的 
22. 1% 和 76% ;草地 和 未 利用 地 向 耕地 转换 的 面积 


分 别 为 2 955 km? 和 3 784 km2 , 占 耕 地 转 和 人 总 面积 
的 42.9% 和 55% ;耕地 和 未 利用 地 向 城市 转换 的 面 
积分 别 为 790 km 和 2 579 km  , 占 耕 地 转 和 人 总 面积 
的 19.3% 和 62. 9% 。2050 年 的 城市 面积 较 2017 年 
将 会 增加 2 倍 , 城 市 的 发 展 将 更 多 地 占用 周边 的 未 
利用 地 和 耕地 ,呈现 为 蔓延 式 扩张 ,各 城市 的 面积 均 
将 有 所 增加 ,增加 最 为 明显 的 是 乌鲁木齐 .吐鲁番 一 
托 克 逊 和 库尔勒 等 地 区 。RCP8.5 情景 高 人 口 和 
GDP 增长 率 , 导 致 耕地 面积 大 幅 增加 ,以 缓和 人 口 
急剧 增长 带 来 的 粮食 生产 压力 ,主要 有 喀什 阿克苏 
和 哈密 等 地 区 。 除 高 海拔 地 区 山顶 终年 积 雪 依旧 存 
在 外 ,其 他 低 海拔 小 范围 的 冰雪 将 会 完全 消失 ,部 分 
冰雪 消融 后 形成 洼地 或 注入 河 湖 , 导致 水 体 总 面积 
增加 。 


#6 2017 年 与 RCP8.5 情景 2050 年 土地 利用 变化 面积 转移 矩阵 


Tab.6 Transition matrix of land use change in 2017 and in 2050 under RCP8.5 scenario /km? 
RCP8.5 情景 2050 年 
上 地 利用 类 型 - 
水 域 林地 草地 耕地 城市 冰雪 未 利用 地 2017 年 合计 
2017 年 水 域 2 842 6 33 4 0 0 11 2 896 
林地 0 914 72 112 2 0 103 1 203 
草地 10 59 188 532 2 955 692 0 1 687 193 935 
耕地 5 2 2 072 46 884 790 0 75 49 828 
城市 0 0 44 0 1 911 0 52 2 007 
冰雪 1 0 27 26 38 2 451 1 659 4 202 
未 利用 地 197 0 7 125 3 784 2 579 137 215 393 229 215 
2050 年 合计 3 055 981 197 905 53 765 6 012 2 588 218 980 483 286 
综 上 所 述 ,3 种 情景 发 展 策略 和 气候 变化 不 同 ， ”为 草地 。 


导致 2050 年 冰雪 面积 明显 减少 ,RCP2.6 、RCP4. 5 
和 RCP8.5 情景 比 2017 年 分 别 减少 29.2% .34.6% 
和 38.4% 。 温 度 上 升 和 城市 发 展 导致 部 分 地 区 水 
体面 积 减 少 ,冰雪 消融 汇 人 部 分 河流 湖泊 或 形成 洼 
地 致使 水 体面 积 增 加 ,总 体 来 说 水 体 总 面积 略 有 增 
加 。 城 市 快速 的 扩张 以 适应 人 口 的 增加 和 经 济 增 
长 ,RCP2.6 和 RCP4.5 情景 下 城市 面积 增加 较 少 ， 
RCP8.5 情景 下 城市 面积 急剧 增加 ,几乎 是 2017 年 
的 3 倍 。RCP2.6 和 RCP4. 5 情景 下 耕地 面积 分 别 
减少 24.3% 和 29.8% ,人 类 的 活动 可 能 占用 更 多 耕 
地 ,使 城市 用 地 增加 ,气候 的 原因 也 使 部 分 耕地 转换 
为 草地 ;RCP8.5 情景 下 耕地 的 增加 7.9% ,以 适应 
城市 以 及 人 口 大 幅度 的 增加 。 前 两 种 情景 林地 面积 
增加 ,而 RCP8.5 情景 下 林地 面积 减少 。 未 利用 土 
地 面积 均 在 减少 ,RCP4.5 情景 减少 最 多 ,为 7.1%; 
未 利用 土地 向 其 他 土地 利用 类 型 转换 ,大 部 分 转换 


3 结论 与 讨论 


本 文 基于 FLUS 模型 ,结合 天 山 山区 53 个 气 
象 站 点 ,31 个 CMIPS 全 球 气候 模式 集合 逐年 平均 
气候 数据 .社会 经 济 和 自然 地 形 数据 对 RCPs 情景 
下 2050 年 天 山 山 区 的 土地 利用 类 型 进行 了 预测 
分 析 : 

(1) 构建 的 FLUS 模型 能 够 较 好 地 模拟 研究 区 
历史 土地 利用 变化 ,具有 较 高 的 模拟 精度 。 

(2) 3 种 情景 下 多 年 平均 气温 升 高 1.67 ~2. 16 
°C ,其 中 春 .夏季 增幅 较 小 , 秋 ,冬季 增幅 较 大 。 天 山 
山区 地 处 西北 干旱 区 ,多 年 平均 降水 量 增加 22 ~ 25 
mm,3 种 情景 下 个 别 年 份 较 基准 期 出 现 降水 减少 现 
象 ,但 总 体 呈 增加 趋势 。 

(3) 与 2017 年 实测 土地 利用 相 比 ,人 类 活动 和 


FOR 区 


气候 变化 导致 2050 年 RCP2. 6 RCP4. 5 和 RCP8.5 
情景 下 的 预 估 结果 相差 各 异 。2050 年 天 山 山 区 冰 
轨 用 地 均 有 减少 且 有 转换 为 草地 或 未 利用 地 之 势 ， 3 
种 情景 下 分 别 减少 了 29.2% 、34. 6% 、38. 4% 。 前 
ee ee 
用 地 和 草地 ,林地 面积 也 有 所 增加 。 而 RCP8. 5 城 
市 用 地 面积 为 2017 年 的 近 3 倍 ,部 分 林地 转换 为 草 
地 。 青 景 水 体 总 面积 略 有 增加 ,未 利用 地 均 呈 
， 31 个 CMIPS 气候 模式 集合 年 均 数 
据 , 没 有 对 各 模式 在 天 山 山区 气候 模拟 的 优 5 
分 析 , 各 模式 模拟 精度 在 气温 和 降水 等 方面 仍 有 差 
异 , 存 在 较 大 的 不 确定 性 。1 km 空间 分 辩 精 度 的 土 
地 利用 图 依旧 会 忽略 地 形 起 伏 植被 变化 和 微 气候 
等 因素 产生 的 差异 ,也 会 带 来 一 定 的 不 确定 性 。 男 
外 , 降 尺 度 方法 的 选择 和 FLUS 模型 参数 的 设置 也 
使 模拟 和 预测 结果 有 一 定 的 误差 。 选 用 其 他 降 尺 度 
的 方法 ,在 实际 多 类 别 土地 利用 变化 模拟 中 ,加 入 其 
他 自然 与 人 类 活动 因素 的 影响 ,将 是 未 来 工作 的 
重点 。 
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Estimation of Future Land Use Change in the Tianshan Mountainous 
Based on FLUS Model 


QIN Qi-rui’”, LI Xue-mei'?, CHEN Qing-wei?, TIAN Ya-lin'? 
(1. Faculty of Geomatics ,Lanzhou Jiaotong University , Lanzhou 730070 , Gansu , China ; 
2. Gansu Provincial Engineering Laboratory for National Geographic State Monitoring , Lanzhou 730070 , Gansu , China; 
3. Transportation Management Office , Jinan Transportation Bureau , Jinan 250014 , Shandong , China) 


Abstract; In this paper, the 31 global climate models were integrated from the Phase 5 of Coupled Model Inter- 
comparison Project (CMIPS) to study the overall prediction of land use change in the Tianshan Mountains under 
three emission scenarios ,i. e. ,the Representative Concentration Pathways (RCPs) 2.6,4.5 and 8.5. In the study, 
the observed climatic data from 53 meteorological stations , socioeconomic data and natural geomorphological condi- 
tions in the Tianshan Mountains were used. The results showed that: T) FLUS model could be used to simulate the 
historical land use change in the study area with high accuracy; @) Compared with the base period from 1970 to 
1999 ,the average annual temperature would be increased by 1.67 -2. 16 °C during 2021 - 2050, and its increase 
would be low in spring and summer but high in autumn and winter. It was revealed that the average annual precipi- 
tation would be increased by 22 -25 mm under RCPs 2.6,4.5 and 8.5. Under the three scenarios , precipitation re- 
duction would occur in some years compared to the base period, but there would be holistically an increase trend; 
(3) Compared with the situation of land use change in 2017 ,it was estimated that, under the RCP2.6,RCP4.5 and 
RCP8.5 scenarios ,the glacier and snow areas in the Tianshan Mountainous in 2050 would be reduced by 29.2% , 
34.6% and 38.4% respectively. The enlargement of urban land area in 2050 would not be significant under 
RCP2.6 and RCP4.5. Some of the cultivated land would be converted into urban land and grassland, and the area 
of woodland would be slightly enlarged. The urban land area under the RCP8.5 scenario would be nearly three times 
of that in 2017. Part of woodland would be converted into grassland. The total area of waters would be slightly en- 
larged ,and the area of unused land will be in a reduction trend. 


Key words: climate change; scenario simulation; land use change; FLUS model; Tianshan Mountainous 


